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Abstract 

We present a model for alloys of compound semiconductors by introducing 

a one-dimensional binary random system where impurities are placed in one 

sublattice while host atoms lie on the other sublattice. The source of disorder 

is the stochastic fluctuation of the impurity energy from site to site. Although 

the system is one-dimensional and random, we demonstrate analytical and 

numerically the existence of extended states in the neighborhood of a given 

resonant energy which match that of the host atoms. 
PACS numbers: 73.20.Jc, 73.61.Ph, 73.20.Dx, 72.20.-i 
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I. INTRODUCTION 



Electron states in random systems have become an active research topic since the gen- 
erality of localization phenomena in one dimension (ID) ffl. Although it is well established 
that almost any nonzero disorder causes exponential localization of all eigenstates in ID sys- 
tems, regardless their energy (see, e. g., Ref. and references therein), there exist several 
exceptions. It is nowadays well known that extended states may arise in random systems 
where disorder exhibits short-range [P|-|T%H or long-range correlations [ID, 211. Spatial corre- 
lation means that random variables are not independent within a given correlation length. 
Suppression of localization by correlations was further put forward for the explanation of 
high conductivity of doped polyaniline as well as transport properties of random semi- 
conductor superlattices ||2T] . 

In this work, we report further progress along the lines in the preceding paragraph. We 
turn ourselves to one of the pioneering works in the field, namely the work of Wu and Phillips 
on polyaniline (see Ref. |7j and references therein). These authors showed that polyaniline 
can be mapped onto a random dimer model that has a set of extended states, originated 
by a resonance of a single dimer defect (two neighbor sites with the same energy). Electron 
states whose energy is close to this resonance turn out to be extended (in the sense that 
their localization length is larger than the system size) when dimers are placed at random 
in the ID system. In this work we show that dimers are not needed to observe extended 
states in ID random systems with short-range correlated disorder. To this end, we built up 
a simple model of semiconductor binary alloy — like ternary III-V compounds — . In these 
alloys (say Al^Ga^a-As), the cation sublattice is occupied by the same atoms (say As) while 
anions (say Al and Ga) are randomly distributed over the other sublattice. We model these 
alloys by considering a ID random binary alloy with two species, referred to as A and B 
atoms hereafter. In order to mimic the disorder present in the anion sublattice, we further 
assume that the site energy of A atoms is randomly distributed from site to site while that 
of B atoms is the same over the entire cation sublattice. As a major point, we demonstrate 
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the occurrence of extended states in the vicinity of the site energy of B atoms in spite of 
the fact that the system is purely ID and random. Thus, we conclude that dimers (or larger 
defects like trimers or n-mers) are not required to observe extended states in ID systems 
with short-range correlated disorder. 

II. MODEL 

We consider a ID binary system where A (B) atoms are placed at odd (even) positions 
of the otherwise regular lattice, whose corresponding site energies are e 2n _i (e 2n ) with n = 
1, 2, ... N, N being the number of unit cells of the alloy. The Schrodinger equation for 
stationary eigenstates ip n (E) is 

{E — e n )ip n + tp n+ i + ip n -i = 0, n = 1,2,..., AT, (1) 

where the E is the corresponding eigenenergy, e n is the site energy and M = 2N is the 
number of atoms in the system. According to our model, site energy at even positions is 
the same and we can set e 2n = without loss of generality. The source of disorder in this 
model arise from the stochastic fluctuations of site energy at even positions. We assume 
that {e 2n -i}^=i is a set of uncorrelated random Gaussian variables with mean value v and 
variance a 2 . Hereafter a will be referred to as degree of disorder. The joint distribution 
function of a realization of disorder is represented by the direct product of single Gaussians. 
Thus 

(e 2n )=v, (4 n ) = v 2 + a 2 . (2) 
The Schrodinger equation (p]) can be written via the 2x2 promotion-matrix P n as follows 
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By iterating this equation we can relate (ipn,ipn+i) and (^on^i) with ip Q = 0: 
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where M n is referred to as the transfer-matrix. We find most convenient to deal with 
the promotion-matrix of the (diatomic) unit cell instead of that corresponding to a single 
atom d^D, namely T n = Pi n Pin-\- F° r rea l E an d e n; the promotion-matrix T n can be 
regarded as an element of the 50(1, 2) group, isomorphic to SL(2, R). It can be cast in the 
following form via the Pauli matrices <r M 



T 

J. , > 



E 



(E - e 2n _!) - 1 



1 2 -^{E - e 2n _i)(73 + -^Tp-^i + * f-B - t -^Y~ ) a 2 , (5) 



where I m denotes the m x m unit matrix. It is easy to demonstrate the following useful 
property T~ x = o^T^o^. The transfer-matrix of the entire system (N unit cells) is obtained 
as M N = IIn=7v ^n- Oseledec's theorem [^2| states that the following limiting matrix T 
exists 



T = lim (mIm n ) 



1/2N 



(6) 



with eigenvalues e 7 . The Lyapunov exponent 7 is nothing but the inverse of the localization 
length A -1 , where A is given in units of the length of the unit cell. 



III. EXISTENCE OF EXTENDED STATES 

In order to find the localization length one should calculate the matrix M^ N M N for large 
N. We will perform this task following the technique developed in Ref. p3[ . By using the 
formula for the decomposition of the product of two spin- 1/2 states into the direct sum of 
scalar and spin-1 states, we have 

{TiYATrX = \{S)%{S)^ + ^K)jAf (<x% (7) 

where 

Af = ^(r/i;rV) (8) 
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is the spin-1 part. Multiplying the expression (|7]) by the left and right by a 2 we have 
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(9) 



Now we should take into account the disorder and calculate the average of V by random 
distribution of e 2n -\ at odd sites 



1 1 N 

(D = ~o 2 ®a 2 + ~(a»a 2 ) <g> {a v a 2 ) I Il( A i) 



(10) 



where Aj is defined by (§). According to Oseledec's theorem p2| , the Lyapunov exponent 
and, correspondingly, the localization length will be given by 



A- 1 = log[£(£)], 



where £(E) is the closest to unity eigenvalue of (Aj), whose elements are 
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- (2 + E A + AEv - 2E 3 v -a 2 -v 2 + {a 2 + v 2 - 4)E' 



-^(E 4 - 2E 3 v + a 2 + v 2 + (a 2 + v 2 - 2)E 2 
—E 3 -v + 2E 2 v - (a 2 + v 2 - 2)E 



- (E 4 + AEv - 2E 3 v - a 2 -v 2 + (a 2 + v 2 - 2)E 2 

- (2 + E 4 - 2E 3 v + a 2 + v 2 + (a 2 + v 2 )E' 
-1 (e 3 + v - 2E 2 v + (a 2 + v 2 )E' 
-2E + E 3 + v - E 2 v 
-i [E 3 - v - E 2 v 
1-2E 2 + 2Ev. 
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Delocalized states have an infinite localization length and, therefore, at some particular 
energy E, the matrix (Aj) should have an eigenvalue equal to one. Hence we obtain the 
following condition for obtaining delocalized states 



det 



Is - (Aj) 



-2a 2 E 2 = 0. 



(13) 



As we see, there is a delocalized state at E = 0. One can calculate the localization length ((TTJ) 
and expand it around E — 0. For v > one gets 

1 2i; 1 / 2 £ 1 / 2 + 0(£ 3 / 2 ), £ > 0. 
Notice that the localization length is asymmetric around the energy of the extended state 
since it scales as ~ E at the left and ~ E l l 2 at the right. The situation is just the opposite 
for v < 0. Remarkably, the prefactor at the left depends on the degree of disorder of the 
alloy but becomes independent of disorder at the right. 

When the degree of disorder vanishes (a = 0), the alloy is simply a diatomic periodic 
chain with site energies and v in each unit cell and, consequently, there are two allowed 
bands. The lower band ranges from v /2 — [(v/2) 2 + 4] 1 ^ 2 up to while the upper band ranges 
from v up to v/2 + [(v /2) 2 + 4]^ 2 for v > 0. Obviously all eigenstates are Bloch functions 
and spread over the entire chain. Localization occurs as soon as a small degree of disorder is 
introduced in the system. But, according to our previous results, the eigenstate with E — 
remains extended. This is clearly seen in Fig. [I], where the inverse of the localization length 
obtained from flTT|) is plotted against E for v = 1 and different degrees of disorder a. In all 
cases the inverse of the localization length is nonzero except at E = 0, where A -1 = 0. This 
suggests the occurrence of a delocalized states at E = 0. 



IV. NUMERICAL RESULTS 



To confirm the above analytical results we have also numerically diagonalized the 
Schrodinger equation (p. We will mainly focus our attention on the normalized density 
of states p(E) and on the degree of localization (inverse participation ratio, IPR) for the 



states at energy E. They are defined respectively as follows |24 
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where the angular brackets indicate an average over an ensemble of disordered linear chains 
and the a^ n is the eigenvector of ([!]) corresponding to the eigenvalue E^ with k — 1, 2, . . . , M . 
Here R is the spectral resolution and 9 is the Heaviside step function. The IPR behaves like 
1 jH for delocalized states spreading uniformly over the entire system on increasing M . In 
particular, the IPR can be exactly computed for the eigenstates of the periodic lattices. In 
doing so we obtain the expected behavior for M — > oo. On the contrary, localized states 
exhibit much higher values. In the extreme case, when the eigestate is localized at a single 
site, the IPR becomes unity. 

We have fixed v — 1 and studied several values of the degree of disorder a, ranging from 
0.1 up to 1.0. The highest value of the degree of disorder considered in the present work 
means that the typical fluctuations of the site energy is of the order of the nearest neighbor 
coupling. The maximum number of atoms in the chain was Af = 1000 (500 unit cells) 
although larger systems were also studied to check that our main results are independent 
of the size. The results comprise the average over 100 realizations of the disorder for each 
given pair of parameters v and a. The spectral resolution was R = 4 x 10 -3 . 

Let us comment the results we have obtained numerically. Figure shows the DOS for 
two different values of the degree of disorder [a = 0.1 and o = 1.0) when the system size 
is M = 1000 and v = 1. The DOS presents the usual U-shape within the bands when the 
degree of disorder is small. The singularities at the edge of the allowed bands are smeared 
out on increasing the degree of disorder except at E = 0, where the divergence remains even 
for the largest degree of disorder (cr = 1.0). This result provides further evidence that the 
states at E = are delocalized. 

The degree of localization (IPR) presents an overall increase when the degree of disorder 
increases, meaning that the larger the degree of disorder, the smaller the localization length. 
This is clearly observed in Fig. ^|, where we show the IPR as a function of the energy for the 
same parameters of Fig. [I]. However, the increase of the IPR strongly depends on the energy, 
being more pronounced close to the center of both allowed bands. Interestingly, the IPR 
at E = becomes independent on the degree of disorder although depends on the system 
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size, as expected. This peculiarity manifest the delocalized character of the states at this 
spectral region. Finally, notice the good correspondence between Fig. [I] (analytical results) 
and Fig. |] (numerical results). 

V. CONCLUSIONS 

In this paper we have considered electron dynamics in a one- dimensional model of binary 
alloy where disorder lies in one of the two sublattices. Although the system is purely one- 
dimensional and random, we have demonstrated analytically the existence of delocalized 
states close to a resonant energy, which match that of the atoms of the other sublattice. 
Numerical results from the evaluation of the DOS and IPR (degree of localization) strongly 
suggest that there exist many states close to the resonant energy that remain unscattered 
in finite systems. 
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FIGURES 

FIG. 1. Inverse of the localization length as a function of energy when v = 1 and a = 0.1 
(dashed line), a = 0.5 (solid line) and a = 1.0 (dotted line). Notice that A — > oo at E = 0. 

FIG. 2. Density of states as a function of energy when the lattice size is M = 1000, v = 1 and 
the degree of disorder is a = 0.1 (dashed line) and a = 1.0 (dotted line). 

FIG. 3. Inverse participation ratio for the same cases shown in Fig. |l]. Notice the overall 
increase on increasing the degree of disorder except at E = 0. 
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